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Abstract 

We demonstrate the newly developed resource for exoplanet researchers - The Ex- 
oplanet Transit Database. This database is designed to be a web application and it 
is open for any exoplanet observer. It came on-line in September 2008. The ETD 
consists of three individual sections. One serves for predictions of the transits, the 
second one for processing and uploading new data from the observers. We use a 
simple analytical model of the transit to calculate the central time of transit, its 
duration and the depth of the transit. These values are then plotted into the ob- 
served - computed diagrams (O-C), that represent the last part of the application. 
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1 Introduction 

Research on extrasolar planets is currently one of the most exciting fields 
in astrophysics. The speculations on the existence of planets orbiting other 
solar-type stars ended fourteen years ago. In 1995 the discovery of the first 
extrasolar planet orbiting a solartype star - the well-known 51 Peg b, was 
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made by iMayor fc Quelozl (119951 ). Since then, the number of known planets 
has been growing quickly. Currently, more than 370 such bodies are knowrPl. 



If a planetary system happens to be oriented in the space so that the or- 
bital plane is close to the line-of-sight to the observer, a planet periodically 
passes in front of the stellar disk. Photometric observation of the transit can 

sical parameters of the planet (e.g., 



Southworth 


2008 or 


Torres et al. 


2008) 



been measured, or that might be measured in the future through precise ob 



servations of transiting planets see IWinnl (120081 ) . At the half of the year 2009 



more than 60 planets with this special orientation were known. 



There are many amateur astronomers all over the world that achieved pho- 
tometric accuracy around the units of percent, which is necessary for quality 
observing of the exoplanet transit. Unfortunately, up to present day, no exo- 
planet light curves' database was available that would accept data from both 
professionals as well as from amateurs. Amateur observers are not constrain 
by telescope scheduling and often have unlimited access to their instrument 
which enables them to gather data over a long period. 



Huge quantity of observations is the key to the search for other planets in 
already known systems. It is importa nt to monitor possible per iodical changes 
in O-C plots of the planets because as lHolman fc Murray! (120051 ) demonstrated 
in their theoretical work, short-term changes of the time of the transit can 
be c aused by the pres e nce of other exo planets or moons in the system, see 
also lAgol et al.l (120051 ). iKipping I (120091 ). On the other hand, potential long- 
term changes in the duratio n of the transit may be the consequence of orbital 
precession of exoplanets as iMiralda-Escudel (120021 ) showed in his theoretical 
work. To perform such effective studies, we need a database which includes all 
available data divided into groups according their quality. 



2 Why ETD 



The Exoplanet Transit Databasd" 2 ! (ETD) came on-line in September 2008 as a 
project of the Variable Star and Exoplanet Section of the Czech Astronomical 
Society. The ETD includes all known transiting planets that have published 
ephemerides. 



1 see the list http://www.exoplanet.eu 

2 http://var2.astro.cz/ETD 
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We have created on-line ETD portal to supply observers with such useful 
information like transit predictions, transit timing variation (TTV), variation 
of depth and duration with availability to draw user observation to the plot. 



Before the ETD, only two other transit databases were available. The Amateur 
Exoplanet Archiv(0(AXA), lead by Bruc e Gary, and the NASA/ IPAC/NExScI 



Star and Exoplanet Database^ (NStED) flvon Braun et all2008f ). AXA strictly 



accepts data only from amateurs. Unfortunately, the quality of light curves is 
very diverse. In spite of this fact, all of the available light curves have the same 
priority grading. On the other hand NStED contains only the light curves that 
have been already published (in the future some amateur light curves from 
AXA should be also accepted into the NStED and the AXA should expire). 

The main goal of the ETD is to gather all available light curves from profes- 
sional and also amateur astronomers (after one year of the existence of the 
database, more than one thousand such records are available). We are search- 
ing for new publications on several open archives to gather all available light 
curves. We also take over data from the NStED, AXA and from our project 
TRESCA0. It is also possible to upload data into the database directly us- 
ing a web-form or it can be added to the database by the administrators. All 
available data are on-line plotted into graphs where we make the provision 
for the quality of the light curve. All graphs (like TTV) can be easily down- 
loaded from the database. It is also possible to download light curves from the 
TRESCA observers and from amateur observers directly from the database. 

While collecting published data to ETD, we accentuate to have fully referenced 
its source. Each transit observation we store full reference with URL pointing 
to the paper or web-page where data were found. When we take over the whole 
light curve we display it only with reference to source of the data there is no 
mention of ETD in the picture. 



3 Parts of the ETD 



The ETD is composed from three sections. The first one - Transit prediction - 
serves for prediction of the transits. The second one - Model fit your data - is a 
web-form for accepting and processing new data. The last section - O-C plots 
- contains the observed - computed diagrams of the central times of transits, 
depths and the transit durations that are generated on-line from the database. 



3 http://brucegary.net/AXA/x.htm 

4 http://nsted.ipac.caltech.edu 

5 http:/ /var2. astro.cz/EN/tresca 
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HD 189733 b 



sefanik observatory Prague 
Htamslau rod dang 




Posi address: 



Statii 



Observing location 
Equipment 
Used filter: 



HD189733 b 



|StanislavPoddany 
|Stefanik observatory. Prague 



| p o d d a ny@ o b servato ry.cz 
|5tefanik observatory. Prague 
ELo-ngitude: 15" | Latitude: 50° 




|Meade LX 2000.40 m + SBIG ST10XME 

ru r B r v cf R r t c am*r 



FHW 10 pixel ejpc^u'E 2C< - 25 s no autoguiding 



Innervation already I 
published: I 

JDmid: |2454705 .30833 +/- 0.0004 
HJDmid: |24E4705.40224 +/- Q.00O4 Jhelcnr = □.□□3911 

Hid transit - UT: |2DDB-G3-26 21:33-:35 

Duration: |l 0-8.2+/- 1.8 ~ minut 



Depth: |Q.0312+/- 0.001 



JO (24547B5 + > 

S - 3.1 mag, ro - 0.93 d 



. . .. ■: « ; 



Fig. 1. A record sample in the ETD after successful processing the light curve. 
3. 1 Transit prediction 



This section of the ETD contains one month prediction of observable tran- 
sits (the starting day is the date two days before current date) and also the 
prediction for the next 365 days for selected exoplanet. Any observer can find 
here the time of the transit start /center/end, duration and the depth of each 
transit for any place in the world. Furthermore the altitude and the cardinal 
point of the object in the sky are displayed for the first contact, mid-transit 
time and the last contact. In the one year prediction window you can also see 
the finding chartGU (15' x 15'). 



3.2 Model-fit your data 



This section describes an user-friendly web-form for uploading and process- 
ing the light curves into the database. To model-fit the transit we assume 
that the observations consist of N relative magnitudes taken at times ti 
(i = 1,2..., N) and the photometry software provided measurement errors cr^ 
computed most likely from Poisson statistics and read-out noise. We model 



6 downloaded from the http:/ /archive.stsci.edu/dss/index.html 
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the dataset by a function 



m(U) =A-2.5\ogF(z[t i ,t ,D,b],p,c 1 ) 
''mean ) + C{U - 

^mean) j 



where F(z,p,ci) is a relative flux from the star due to the transiting planet. 
We assume that the planet and the star are dark and limb darkened disks, 
respectively, with radius ratio p = Rp / R* and that the planet is much smaller 
than the star, p < 0.2. The projected relative separation of the planet from 
the star is z. Limb darkening of the star is modeled by the linear law with the 



coefficient c\. We employ the occultsmall routine of iMandel fc Agoll (120021 ) as 
our F(z,p, ci). We checked that the small planet approximation, p < 0.2, does 
not produce significant differences from the full model (at least for the typical 
values of p and having in the mind the typical quality of the photometry) 
and is much faster to compute, which is the most important factor for on-line 
processing. 

We model the planet trajectory as a straight line over the stellar disk with 
impact parameter b = a cos I/R*, where a is a semi-major axis and / is the 
orbit inclination. The mid-transit occurs at to and the whole transit lasts D. 
Based on these assumptions we can compute z[ti, t , D, b] for every ij. 

Variable A in the equation ([I]) descibes the zero-point shift of the magni- 
tudes, while B and C describe systematic trends in the data. Linear and 
quadratic terms are computed with respect to the mean time of observations 
Wan = J2 U/N to suppress numeric errors. We do not employ any explicit cor- 
rection for air-mass curvature as we think a generic second-degree polynomial 
is sufficient in most cases. 

We used the Leyenberg -Marquardt non-linear least squares fitting algorithm 



from lPress et al.l (119921 ). procedure mrqmin. The algorithm requires initial val- 
ues of parameters and partial derivatives of the fitted function. We take the 
initial values from literature (except for c±, see below). We compute all par- 
tial derivatives of the equation ([1]) analytically, except for dFj dz, dF/ dp and 
dF/dci which were comp uted numerically using Ridders' method (procedure 



dfridr of iPress et al.lll992f ) 



The search for the optimal parameters is done by iterating the fitting proce- 
dure until the A% 2 (between fits) does not change significantly. Usually, with 
good initial values, about ten iterations are sufficient. Then the error bars 
Oi are re-scaled to make the final x 2 = N — g, where g is a number of free 
parameters, and we re-run the fitting procedure to obtain final errors of the 
parameters. Original photometric errors are usually underestimated and this 
procedure yields more reasonable errors of the output parameters. 
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Table 1 

The comparison of our results with IWinn et al.l (|2007l ) . 



HD189733b 


ETD 


Winn et al. C2007) 


central time [HJD] 


2453988.80333 (12) 


24543988.80331 (27) 


duration [min.] 


106.01 ± 0.50 


109.62 ± 1.74 



In the optimal circumstances, one would consider all variables in equation ([I]), 
namely A, B, C, to, D, b, p and c% free parameters. However, these parameters 
are correlated to some extent and noisy photometry from a small amateur 
telescopes does not permit recovery of all of them. We need to fit the zero- 
point shift A and in most cases also a linear systematic trend B. Our primary 
goal is to get the central time of the transit to, duration D and depth of the 
transit 5. Hence, we set to and D as free parameters, by default. However, for 
a limb-darkened star, the depth of the transit S is determined by radius ratio 
p, impact factor b and limb darkening coefficient c±. Primarily, the depth of 
transit 5 is governed by the radius ratio p and we set it as a free parameter. 
The parameters b and c\ affect the depth and the shape of the transit to a 
lesser extent and from noisy amateur data we couldn't retrieve meaningful 
values for the two parameters simultaneously with p. Therefore, we hold b and 
C\ fixed. We either compute b from orbital parameters of the planet and from 
the radius of the star or take the value from the literature. The situation is 
more complicated for limb darkening because c\ should be different for every 
photometric filter. We decided to keep c\ fixed at a rather arbitrary value 
Ci = 0.5 in all cases. We experimented with values from 0.2 to 0.9 and found 
that the effect on other parameters is rather negligible, usually smaller than 
the error bars. The export value of the depth is then evaluated as 



-2.5 log 



minF (z,p, C\) 



(2) 



At Fig. [T] you can see the example of the record in the ETD after successful 
processing of the light curve. 

We tested our al gorithm using the data of the exoplanet HD189733b that 



were published by lWinn et al.l (120071 ). When we fitted this precise data using 
our code we obtained the value of the mid-transit time which was in excellent 
agreement with the published results. The duration of the transit wasn't inside 
the error bars. The errors of our fit are lower than those in the paper where 
the MCMC simulations were used (Tab. [I]). We think that this excess is due 
to red noise which we did not take into account in our calculation and because 
of impact parameter that is fixed in our case. We also made a test with the 
data from AXA and we obtained similar results to the results presented in 
their archive (Tab. [2]). 



6 



Table 2 

The comparison of our results with the AXA database. 



HD189733b 



ETD 



AXA 



central time [HJD] 


2454705.40228 (41) 


2454705.4023 (5) 


duration [min.] 


102.09 ± 1.60 


98.4 ± 1.8 


depth [mag.] 


0.0287 ± 0.0006 


0.02895 ± 0.00080 


TrES-l b 


Exoplanet Transit Database 


O-C vs EPOCH 




-170 -90 -10 70 150 

EPOCH (yeers 2003 - 2008, 45 records) 



Fig. 2. The central transit time O-C plot for the exoplanet TrES-l b in the ETD. 
3.3 O-C plots 



The O-C plots section contains the diagrams of the central transit time, the 
duration of the transit and its depth as a function of time. All available data 
for the selected planet including the error bars are visible in these graphs 
(Fig. [2]). The quality of individual light curve is indicated by the size of the 
dot. The records in the ETD are divided into 5 groups according to their data 
quality index DQ. While computing DQ index of the light curve, the following 
relation is used 



8 r 



(3) 



where a is a temporary data quality index, S is the mean absolute deviation 
of the data from our fit and p = N/l means the data sampling, where I is the 
length of observing run in minutes. 
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Table 3 

The distribution of the quality of the light curves according to their DQ index 
(example of the light curves see Figure [3]) . 

DQ index 1 2 3 4 5 

threshold a > 9.5 9.5 > a > 6.0 6.0 > a > 2.5 2.5 > a > 1.3 1.3 > a 
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Fig. 3. Examples of the light curve with the different DQ index. 

The number a is further transformed for better lucidity to the scale from 1 
to 5 where 1 presents the best quality data and the value 5 the worst data 
(Tab. [31 Fig. [3]). These thresholds are used only if whole transit is observed or 
when we take over whole light curve (AXA, TRESCA). If some part (egress or 
ingress) is missing in the dataset, the observation automatically gets the worst 
DQ index and in the notice column in the summary table the notice " Only 
partial transit" is generated. When we take over only the results of midtransit 
time, transit depth and length of the transit (not the whole light curve) we 
usually give the DQ equal 1. 



3.3.1 How to download data 

In the O-C plots section of the ETD you can also download the data for 
your further studies. If there is displayed whole light curve you can download 
it using link in the DQ column. Whole table including the O-C residuals 
can be downloaded directly using link "Show data as ASCII table separated 
by semicolon" below the table. If the transit observation that you used in 
your next study was captured to ETD from literature, you can find reference 
in column "Author & REFERENCE", so the source paper should be cited 
in common way. If transit observation was published in some on-line source 
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(AXA, TRESCA), ask observers for a permission and other useful comments 
about the data. We can supply you with e-mail contact to observer if necessary. 



4 Future development and discussion 

We have created the on-line portal to supply observers with information like 
transit predictions, TTV, variation of depth and duration plots with avail- 
ability to draw user observation to the plot. We think that main part of the 
database is now ready to use and should be useful tool for community. 

In future we plan some improvements of the database to be more user-friendly. 
We also plan to implement limb darkening tables into our fits. Further in the 
future we would like to develop a (semi) automatic procedure for searching of 
the transit timing variations in the O-C diagrams. 
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